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The angle of repose for the flow of granular materials 
in a half-filled rotating drum is studied by means of exper- 
iments and computer simulations. Particles of different ma- 
terial properties are used to investigate the effects of the end 
caps on the angle of repose. By fitting the numerical results 
to an exponentially decaying function, we are able to calcu- 
late the characteristic range, £, of the influence of the wall. 
We found that £ scales with the drum radius but does not 
depend on either the density or the gravitational constant. 
For increasing particle diameter, finite size effects are visible. 
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I. INTRODUCTION 

The behavior of granular materials is of great techno- 
logical interest [Q and its investigation has a history of 
more than two hundred years. When granular materials 
are put in a rotating drum, avalanches are observed along 
the surface of the granular bulk In industrial pro- 

cesses, such devices are mostly used for mixing different 
kinds of particles. However, it is also well known that 
particles of different sizes tend to segregate in the radial 
and axial directions [^-pd| . 

Recently, the particle dynamics of granular materials 
in a rotating drum has been described by using quasi 
two-dimensional systems, tracking individual grains via 
cameras and computer programs Extensive numeri- 
cal studies have also reproduced and predicted many of 
the experimental findings |15-21|. The segregation and 
mixing process depends on many parameters, such as 



size |p|, pi|, s hape |19[, mass jlfj, frictional forces, angular 
velocity] 2lJ , filling level of the drum p3), etc. The angle 
of repose of the material also depends on the parameters 
and it was argued that either the dynamic or static angle 
difference of the materials in the drum influence the axial 
segregation process [[||] |n],|llj] . 

In this article, we investigate experimentally the de- 
pendence of the dynamic angle of repose on the rotation 
speed of a half-filled drum for particles of different ma- 
terial properties. It is found that the angle is up to 5 
degrees higher at the end caps of the drum due to bound- 
ary friction. Using a three-dimensional discrete element 
code, we are able to quantify this boundary effect and 
discuss its dependence on gravity, particle size and den- 
sity. 



II. EXPERIMENTAL RESULTS 

An acrylic cylinder of diameter 6.9 cm and length 49 
cm was placed horizontally on two sets of roller supports 
and was rotated by a well-regulated electronic motor. 
The material used was mustard seeds which are relatively 
round of average diameter about 2.5 mm, and have a co- 
efficient of restitution, estimated from a set of impact 
experiments, of about 0.75 |23|. A set of experiments 
were conducted to measure the angle of repose in dif- 
ferent flow regimes. For a small rotation speed, f2, in- 
termittent flow led to a different angle before and after 
each avalanche occurred, called the starting (maximum) 
and stopping (minimum) angle, respectively. For a larger 
rotation speed these intermittent avalanches became a 
continuous fiat surface and thus enabled to define one 
angle of repose defined as the dynamic angle of repose as 
shown in Fig. [j]a. When fl increases, the fiat surface de- 
forms with increasing rotation speeds and develops a so- 
called S-shape surface for higher rotation speeds, shown 
in Fig. [l]c. The deformation mostly starts from the lower 
boundary inwards and can be well approximated by two 
straight lines with different slopes close to this transi- 
tion, sketched in Fig. [l]b. For all measurements in this 
regime, we took the slope of the line to the right which 
corresponds to the line with the higher slope. 

The average maximum and minimum angles of repose 
for the intermittent avalanches were found to be about 36 
and 32 degrees, respectively, see Fig. ||. There seems to 
be a rather sharp transition from intermittent to continu- 
ous avalanches, which happens around = 4 rpm. For 
greater than 4 rpm where the avalanches are continuous, 
the mustard seed data indicate a linear dependence of 
the dynamic angle of repose on the rotation speed which 
differs from the quadratic dependence found by Rajchen- 
bach §. 

We also investigated the dynamic angle of repose for 
different particle diameters and materials in the contin- 
uous regime in more detail using a 27 cm long acrylic 
cylinder of diameter 2R = 6.9 cm. For a given rota- 
tion speed, f2, the dynamic angle of repose was measured 
four times at one of the acrylic end caps and the aver- 
age value with an error bar corresponding to a confidence 
interval of 2a, where a is the standard deviation of the 
data points, was then calculated. First we used mustard 
seeds of two different diameters, namely 1.7 mm (black) 
and 2.5 mm (yellow), with a density of 1.3 g/cm 3 . The 
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latter were the same that were used to produce Fig. |2j. 
We varied the rotation speed, fi, from 5 rpm to 40 rpm 
and took the higher angle in the S-shaped regime which 
exists for higher rotation rates, see Fig. 0b. Both data 
sets are shown in Fig. |^ for black (•) and yellow (o) seeds. 
The figure also illustrates the transition to the S-shaped 
regime which occurs at the change of slope, e.g. at around 
11 rpm for the smaller seeds and around 16 rpm for the 
larger seeds. One also notes that the dynamic angle of 
repose is much higher for the larger particles in the low 
frequency regime. For values of f2 > 15 rpm in the S- 
shaped regime, the difference in the dynamic angle of 
repose for the two different types of mustard seeds de- 
creases with increasing f2, and both curves cross around 
30 rpm giving a slightly higher angle for the smaller seeds 
with the highest rotation speeds studied. 

We applied the same measurements to two sets of glass 
beads having a density of 2.6 g/cm 3 . The smaller beads 
had a diameter of 1.5 mm with no measurable size dis- 
tribution, whereas the larger beads had a diameter range 
of 3.0 ± 0.2 mm. Both data sets are shown in Fig. |] 
for small (•) and large (o) beads. It can be seen from 
this figure that the transition to the S-shaped regime oc- 
curs at around 16 rpm for the smaller beads and around 
24 rpm for the larger beads. In general, we found that 
the small particles exhibit the S-shaped surface at lower 
values of f2 than the large particles. The angles of re- 
pose are, in general, lower for the glass beads compared 
to the mustard seeds which we attribute to the fact that 
the mustard seeds are not as round as the glass beads 
and rotations of the mustard seeds are therefore more 
suppressed. The coefficient of friction is also higher for 
mustard seeds. 

There are two striking differences when comparing 
Figs. |^ (glass spheres) and || (mustard seeds). For ro- 
tation speeds, f2, lower than 15 rpm, the small and large 
glass beads have the same dynamic angle of repose which 
agrees with the findings in p4| , whereas the dynamic an- 
gle of repose is significantly higher (3 to 4 degrees) for 
the larger mustard seeds compared to the smaller ones. 
For rotation speeds, Q, higher than 15 rpm, the smaller 
glass beads show a higher dynamic angle of repose than 
the larger glass beads, and this angle difference increases 
with increasing rotation speed. For mustard seeds, the 
difference in the dynamic angle of repose between the 
smaller and the larger particles decreases with increasing 
il and the smaller seeds only show a higher angle for the 
highest rotation speeds studied. Both Fig. || and Fig. |] 
seem to indicate that the increase in the dynamic angle 
of repose with rotation speed, f2, in the S-shaped regime 
is larger for the smaller particles. 

All the above angle of repose were measured by looking 
through one of the acrylic end caps. In order to study 
the boundary effect of these end caps on the dynamic an- 
gle of repose, we performed Magnetic Resonance Imag- 
ing (MRI) measurements. This technique of studying 
non-invasively the flow properties of granular materials 
was first used by Nakagawa et al. p3| and is in addi- 



tion explained in more detail in ref. M. We used the 
large mustard seeds, which had an average diameter of 
2.5 mm. The dynamic angle of repose was measured 
based on the concentration data which was averaged in a 
thin cross-sectional slice in the middle of the cylinder far 
away from the end caps. It is shown in Fig. ^| as function 
of the rotation speed by the open circles (o). We re- 
stricted the measurement to the flat surface regime, and 
all data points then lie approximately on a straight line. 
On the other hand, non-MRI data were measured at the 
end caps and are shown as stars (*) in Fig. || The con- 
sistently higher dynamic angle of repose at the end cap 
indicates the significance of the friction between particles 
and the boundary wall. We also found that the S-shape 
regime seems to start earlier at the end caps due to the 
additional wall friction. 



III. SIMULATION TECHNIQUE 

When measuring the dynamic angle of repose for dif- 
ferent materials at the end caps, we found that the angle 
is always lower in the middle of the drum and the influ- 
ence of the end caps seems to be rather short range, the 
angle drops to the value in the middle of the drum within 
a few centimeters. We are particularly interested in the 
dependence of this length scale on the particle diameter 
and density and on the gravitational force. Since we only 
have a limited number of particle diameters and density 
available in the experiment, we use three-dimensional dis- 
crete element methods, also known as granular dynamics, 
to overcome this problem. 

Each particle i is approximated by a sphere with radius 
r,;. Only contact forces during collisions are considered 
and the particles are not allowed to rotate. Since the 
mustard seeds are slightly aspherical, they rotate much 
less than glass beads. This was the motivation for our 
non-rotating assumption. The forces acting on particle i 
during a collision with particle j are 

Ffj = -Y (n + rj - fyn) - 7n%n (1) 

in the normal direction (n) and 

F^ = -m^( ls v ir s{t)^\Fl)\) • (2) 

in the tangential direction (s) of shearing. In Eqs. ([!]) 
and (H), 7n and 7 S represent a dynamic friction force 
in the normal and tangential direction, respectively, 
represents the vector joining both centers of mass, Vij 
represents the relative motion of the two particles, and 
Y is related to the Young's Modulus of the investigated 
material. Dynamic friction in the model is defined to be 
proportional to the relative velocity of the particles in 
the tangential direction. 

During particle-wall contacts, the wall is treated as 
a particle with infinite mass and radius. In the normal 
direction, Eq. (|lj) is applied, whereas in the tangential 
direction, the static friction force 
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Ffj = -min(fc s J Vij ■ 8(t)dt,n\FS\) 



(3) 



is used. This was motivated by the observation that when 
particles flow along the free surface, they dissipate most 
of their energy in collisions and can come to rest in voids 
left by other particles. This is not possible at the flat 
drum boundary. In order to avoid additional artificial 
particles at the walls which would make the simulations 
of three-dimensional systems nearly infeasible, we rather 
use a static friction law to avoid slipping and allowing 
for a static surface angle when the rotation is stopped. 
Both tangential forces are limited by the Coulomb cri- 
terion which states that the magnitude of the tangential 
force cannot exceed the magnitude of the normal force 
multiplied by the friction coefficient p. For particle- 
particle collisions we use p = 0.2, and for particle-wall 
collisions, p, w = 0.4. In order to save computer time, we 
set Y to 6 • 10 4 Ba m which is about one order of mag- 
nitude softer than vulcanite but the maximal overlap of 
two particles is at most 0.3% of the sum of their radii, 
which is still realistic. This gives a contact time during 
collisions of 1.1 ■ 10~ 4 s. The coefficient of restitution 
for wall collisions is set to 0.77 which is within the error 
bar of the experimentally measured value of 0.75, see sec- 
tion 0. In experiments with spherical liquid-filled parti- 
cles, we found only a weak dependence of the restitution 
coefficient on particle size and therefore used a normal 
force law, Eq. ([!]), that would make the restitution coef- 
ficient independent of particle size. When the same type 
of force law is applied to particle-particle collisions, it 
gives a normal restitution coefficient of 0.56. A discus- 
sion of the different force laws is given in ref. |25| and 
a review of different applications using granular dynam- 
ics is given in |p6| . Even though detailed experiments for 
binary collisions of particles were performed, the force re- 
lations before and after a collision depend on the material 
and the aspherity of the particles |^7j and since these two 
quantities were not available for mustard seeds, we can 
only take the published values as a guideline. The nu- 
merical parameters were fine adjusted by comparing the 
experimentally determined dynamic angle of repose for 
2.5 mm mustard seeds (p = 1.3 g/cm 3 ) with the simu- 
lation results over the fi-range from 8 to 35 rpm. The 
radius of the drum was chosen as R = 3.5 cm. Both data 
sets are shown in Fig. ^. Also shown as a solid line in 
Fig. ^| is the theoretical result based on a model by Zik 
et al. |p4| . They started from the equilibrium condition 
for the surface flow j in a laminar and thin layer inclined 
with an angle O: 
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where p denotes the particle density, g gravity, rj the 
constant viscosity and Go = arctan p. The cut-off 
depth, ho, corresponds to a constant pressure value of 
Po = h gcosQ. A second expression for the surface flow 
in a half-filled drum can be obtained by looking at mass 
conservation, pi: 
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(5) 



where r measures distance from the drum center along 
the free surface. Equating expressions (|J) and (|J) and 
using the relation tanO = = y', where y(x) measures 
the height of the top surface particle along the surface 
and (cosG) -1 = yl + (y 1 ) 2 , we obtain 

(y') 3 ~ {y'? tan Go + y' + cft(y 2 + x 2 - R 2 ) = tan9 

(6) 

with c = I^t". Corrections to this model were recently 

proposed by Khakhar et al. |2Sj], but they lead to the 
same equations for the dynamic angle of repose as above 
in the case of a half-filled drum due to the symmetry of 
the thickness of the fluidized layer for shear flow. Solving 
for y' at the origin (drum center), the only one of the 
three roots with no imaginary part reads 



y = A + 



1/3 -n 



(7) 



where 3/2 = tanOo, 



B = p(l + p 2 ) + -dlR 2 and 
D = 3(/i 2 + i) 2 + p(l + p 2 )cflR 2 + (^cilR 2 ) 2 . 

We integrated Eq. (j6|) numerically and checked that the 
theoretical profile has a similar shape for different ro- 
tation speeds as the numerical data. We adjusted the 
parameter c using the experimental data points and the 
best fit was obtained for c = 0.0111 ± 0.0001 s/m 2 and 
p = 0.51 => Co = 27°. It is remarkable how well the 
theoretical curve fits the data points from the experi- 
ments and the numerical simulations, as shown in Fig. |^. 
The theoretical curve is very close to the arcus-tangent 
curve proposed by Hager et al. |^] and also found in 
two-dimensional simulations over a wide ^-regime [ [30"| . 
For low rotation speeds, Q < 8 rpm, the experiment is 
very near the discrete avalanche regime, and therefore 
the simulations where we have used only dynamic fric- 
tion for the particle-particle interactions and the theory 
where a steady flow is assumed tend to deviate slightly 
from the experiment. Rajchenbach experimentally found 
the relation ~ |C — C c | m with m = 0.5 leading to 
C = C c + afl 2 J3| which gives an increasing slope for in- 
creasing in the graph, whereas the experimental data 
points in our Fig. ^ suggest a decreasing slope with in- 
creasing tt. To illustrate this point further, we replot- 
ted all of our available experimental data points of the 
large mustard seeds measured at the end caps, taken from 
Fig. H above, in the same fashion as Rajchenbach and ob- 
tained a scaling exponent of m=0.87 using C c = 34.1°. 
This has to be compared with m=0.5 found by Rajchen- 
bach and m=0.7 given by the numerical prediction by 
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Tang and Bak pl||. The SI range of Rajchenbach is 
smaller than the one investigated by us and we speculate 
that his finding is valid close to the transition point to 
the continuous flow regime where the quadratic fit works 
rather well. We included this in Fig. || as dotted line 
using a best quadratic fit for the value of S7 < 20 rpm. 



IV. BOUNDARY EFFECT ON SURFACE ANGLE 

As shown in Fig. g the dynamic angle of repose 
of granular material in a rotating drum is significantly 
higher at the end caps than in the middle. For the ex- 
perimentally investigated particle sizes of the order of 
millimeter, the effect was visible up to a few centimeters. 
But this length scale might depend on the particle di- 
ameter and density and also on external parameters like 
gravity. 

Using the above described technique, we simulated ex- 
tended three-dimensional, half- filled drums. For a parti- 
cle size of 2.5 mm and a rotation speed of 20 rpm, the 
time-averaged angle, denoted by (9(z)), as function of 
the position along the rotation axis is shown in Fig. |t]. 
Each point with corresponding error bar stands for a 
weighted average over the nearest neighbours. In order to 
study the characteristic length, £, of the boundary effect, 
we fit all data points by the relation 



(6(z)) = e oc + AQ (V^ + e-( L - z ^) 



(8) 



where L stands for the length of the drum, 9^ for the 
dynamic angle of repose far away from the boundaries 
and AO for the angle difference between the value at the 
boundary and 0oo. For the curve shown in Fig. |?], the 
corresponding values are AO « 4°, which is the same 
value given in Fig. || for the MRI experiment, £ = 3.19 ± 
0.25, L = 20.6 cm and 9^ = 35° ± 0.2°. We tested 
Eq. (^) against the simulation results for different drum 
lengths of L/2, L/4 and L/8 and found a remarkably 
good agreement. In the last two cases, the value for 9^ 
is never reached in the middle of the drum due to the 
boundary effects. 

Using drums that are at least 64 particle diameters 
long, we studied the dependence of 9oo (middle) and the 
angle at the end caps on the particle diameter d at the 
same rotation speed of 20 rpm. The simulation values in 
the tangential direction were chosen in such a way that 
the normalized tangential velocities before and after im- 
pact were independent of the particle diameter d. The 
ration Y/k s was set to 3.1, a value which gives for ac- 
etate spheres a good agreement of simulations pa] and 
experiments |27]] . The value of 7 S was chosen sufficiently 
high to give a similar behavior in particle-particle- and 
particle-wall-collisions in the sliding regime. The results 
for the surface angle along the rotation axis are shown in 
Fig. ||, which illustrates that the angle increases with in- 
creasing particle size in agreement with the mustard seed 



experimental results given in Fig. |^. The angle difference 
of around 4 degrees, which seems to be independent of 
the particle size, also agrees with the experimental find- 
ings. In other experiments, different dependencies were 
observed: das Gupta et al. (6) mostly found a higher an- 
gle for smaller particles using sand grains and Hill and 
Kakalios |l(| measured higher angles for smaller parti- 
cles when using sand and glass particles although it was 
also possible to get no angle difference for certain size 
ratios of glass spheres. The latter was also found by Zik 
et al. |24[ ], whereas Cantelaube |32| did not find a clear 
trend when using discs in a quasi two-dimensional drum. 
What causes the different behaviour is not clear at the 
moment and a more detailed analysis would be desirable 
but is beyond the scope of this article. It is necessary to 
use appropriate values for the simulation parameters to 
quantitatively model a desired system, which is why we 
gathered as much information for the mustard seeds as 
possible. A arcus-tangent fit which gave a smaller mean 
deviation than a parabolic fit was added to Fig. || to 
guide the eye. Changing the density, p, of the particles 
or the gravitational constant, g, has a dramatic effect on 
the angle of repose: for the latter quantity this is shown 
in table | and a similar behaviour was seen in recent ex- 
periments p3| . For both p and <?, an increase in value 
corresponds to an angle decrease. When a hydrostatic 
pressure, pq ~ g, is assumed, the data for g < 30 m/s 2 
can be well described by Eq. (0). The lower g becomes 
the more pronounced the S-shaped surface becomes, and 
in the limit g — > 0, the transition to the centrifugal regime 
takes place at 



y/2i 



R sin 9o 



where 9o denotes the average angle in the limit SI — > 
[15 3(J. In our case, for SI = 20 rpm the transition to 
the centrifugal regime occurs at g 0.45 m/s 2 , which is 
in perfect agreement with the numerical findings. Even 
though we studied more than one (two) orders of mag- 
nitude in p (g), we could not obtain an accurate infinite 
value limit. 



V. RANGE OF BOUNDARY EFFECT 

In order to study the range of the boundary effect, we 
extracted figures similar to Fig. [?] from our simulations 
and varied the drum length, L, and radius, R, the particle 
diameter, d, and density, p, the gravitational constant, g, 
and the rotation speed, SI. The data points for the dy- 
namic angle of repose as a function of position along the 
rotation axis were fitted by Eq. (||) giving the charac- 
teristic length, (, of this run. As expected, £ did not 
vary when the length of the drum or the rotation speed 
was changed, but surprisingly, the characteristic length, 
C, did not change when the density of the particles or the 
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gravitational constant was changed by more than one 
order of magnitude, even though the dynamic angle of 
repose strongly depends on both as shown in table | for 
the latter quantity. 

Based on the definition of £ in Eq. (||), one might spec- 
ulate that C ~ R since the gradient of the slope along the 
rotational axis of the surface should be a material prop- 
erty, i.e. it should not depend on the geometry. The angle 
of repose is independent of the drum radius, R, and there- 
fore the height difference between the surface at the end 
cap and the surface in the middle of the drum must be 
proportional to R. This leads to £ ~ R which is indeed 
the case, and we show in Fig. |9| the dimensionless charac- 
teristic length, C/R, as function of dimensionless particle 
diameter, d/R, for three different drum radii. Below a 
critical diameter, d c , £ seems to be independent of the 
particle size, and we propose the following relation 

r { aR , if d < d c , Q , 

C ~ \ aR + p(d-d c ) , if d>d c [) 

where, in our case, a — 0.28 and (3 = 3.13. The critical 
particle diameter d c « 0.14 7? and it seems to decrease 
slightly with increasing drum radius. Therefore, particles 
in the fluidized zone with d < d c might be describable by 
a continuum model. For particles with d > d c , we have 
to take finite size effects into account. 

For comparison, we replot in Fig. [To] the data from 
Fig. H by showing the characteristic length, £, made di- 
mensionless by the average particle diameter, d. A strong 
decrease in £/d and a later saturation is clearly visible for 
increasing particle diameter. The solid line is an expo- 
nential fit which matches all data points very well but 
it can only serve as a guide to the eye since it does not 
reproduce the right value in the limit d — * 0. 

VI. CONCLUSIONS 

We have investigated the dynamic angle of repose, O, 
in a three-dimensional rotating drum in the continuous 
flow regime. By choosing different materials and particle 
diameters, we discussed the ^-dependence of 9 for glass 
beads and mustard seeds of two different sizes. In the low 
rotation speed regime, both types of glass beads showed 
the same angle of repose, whereas the angle was higher 
for the larger mustard seeds. Using MM techniques, we 
could quantify, for the large mustard seeds, the angle dif- 
ference between the middle and the end of the drum and 
its fi-dependence. In all cases, either a linear or an arcus- 
tangent dependence of on Q was found. In order to 
investigate the range, £, of the boundary effects, we used 
a three-dimensional discrete element code and fitted the 
averaged angle along the rotation axis to two exponen- 
tially decaying functions. We found that ( scales linearly 
with the drum radius. On the other hand, it does not 
depend either on the particle density or the gravitational 



constant, even though the surface angle changes drasti- 
cally with these quantities, or on the rotation speed of the 
drum. A detailed analysis of the dependence of the char- 
acteristic length, C on the particle diameter, d, revealed 
that £ is independent of d for small particle diameters 
but shows finite size effects for larger d. 

ACKNOWLEDGMENTS 

CMD and GHR gratefully acknowledge financial sup- 
port by the Deutsche Forschungsgemeinschaft and par- 
tial funding for their stay in Golden by the Particulate 
Science & Technology Group of the Colorado School of 
Mines. JLM and MN are supported in part by NASA 
through contract number NAG3-1970. We also would 
like to thank Susan McCaffery for a critical reading of 
the manuscript. A generous grant of computer time on 
the Cray T3E at the Forschungszentrum Jiilich made the 
numerical simulations possible. 



[1] H. M. Jaeger, S. R. Nagel, and R. P. Behringer, Physics 

Today 4, 32 (1996). 
[2] P. Evesque and J. Rajchenbach, C.R. Acad. Sci. Paris, 

Serie II 307, 223 (1988). 
[3] J. Rajchenbach, Phys. Rev. Lett. 65, 2221 (1990). 
[4] M. B. Donald and B. Roseman, Br. Chemical Engineering 

7, 749 (1962). 
[5] J. Bridgewater, Powder Technol. 15, 215 (1976). 
[6] S. Das Gupta, D. V. Khakar, and S. K. Bathia, Chem. 

Engineering Science 46, 1513 (1991). 
[7] M. Nakagawa, Chem. Engineering Science-Shorter Com- 
munications 49, 2544 (1994). 
[8] E. Clement, J. Rajchenbach, and J. Duran, Europhys. 

Lett. 30, 7 (1995). 
[9] F. Cantelaube and D. Bideau, Europhys. Lett. 30, 133 

(1995). 

[10] K. M. Hill and J. Kakalios, Phys. Rev. E 49, 3610 (1994). 
[11] K. M. Hill and J. Kakalios, Phys. Rev. E 52, 4393 (1995). 
[12] K. M. Hill, A. Caprihan, and J. Kakalios, Phys. Rev. 

Lett. 78, 50 (1997). 
[13] M. Nakagawa, S. A. Altobelli, A. Caprihan, and E. 

Fukushima, Chem. Engineering Science 52, in press 

(1997). 

[14] K. M. Hill, A. Caprihan, and J. Kakalios, Axial Seg- 
regation of Granular Media Rotated in a Drum Mixer: 
Pattern Evolution, preprint, 1997. 

[15] O. R. Walton and R. L. Braun, in Proceedings of the 
joint DOE/NFS workshop on FLOW OF PARTICU- 
LATES AND FLUIDS, edited by M. Roco and S. Plasyn- 
ski (Ithaca, New York, 1993). 

[16] G. H. Ristow, Europhys. Lett. 28, 97 (1994). 

[17] G. Baumann, I. Janosi, and D. E. Wolf, Europhys. Lett. 
27, 203 (1994). 



5 



[18] G. Baumann, I. M. Janosi, and D. E. Wolf, Phys. Rev. 

E 51, 1879 (1995). 
[19] V. Buchholtz, T. Poschel, and H. J. Tillemans, Physica 

A 216, 199 (1995). 
[20] G. H. Ristow, Europhys. Lett. 34, 263 (1996). 
[21] C. M. Dury and G. H. Ristow, J. Phys. I France 7, 737 

(1997). 

[22] G. Metcalfe, T. Shinbrot, J. J. McCarthy, and J. M. Ot- 

tino, Nature 374, 39 (1995). 
[23] M. Nakagawa, S. Altobelli, A. Caprihan, E. Fukushima, 

and E.-K. Jeong, Experiments in Fluids 16, 54 (1993). 
[24] O. Zik, D. Levine, S. G. Lipson, S. Shtrikman, and J. 

Stavans, Phys. Rev. Lett. 73, 644 (1994). 
[25] J. Schafer, S. Dippel, and D. E. Wolf, J. Phys. I France 

6, 5 (1996). 

[26] G. H. Ristow, in Annual Reviews of Computational 
Physics I, edited by D. Stauffer (World Scientific, Sin- 
gapore, 1994), p. 275. 

[27] S. F. Foerster, M. Y. Louge, H. Chang, and K. Allia, 
Phys. Fluids 6, 1108 (1994). 

[28] D. V. Khakhar, J. J. McCarthy, T. Shinbrot, and J. M. 
Ottino, Phys. Fluids 9, 31 (1997). 

[29] W. Hager, S. Linz, and P. Hanggi, in Friction, Arching, 
Contact Dynamics, edited by D. E. Wolf and P. Grass- 
berger (World Scientific, Singapore, 1997), p. 317. 

[30] C. M. Dury, G. H. Ristow, and M. Nakagawa, in Powders 
& Grains 97, edited by R. P. Behringer and J. T. Jenkins 
(Balkcma, Rotterdam, 1997), p. 499. 

[31] C. Tang and P. Bak, Phys. Rev. Lett. 60, 2347 (1988). 

[32] F. Cantelaube, Ph.D. thesis, University of Rennes I, 
Rennes, France, 1995. 

[33] R. Peralta-Fabi, E. Morales, and V. Romero-Rochin, in 
Powders & Grains 97, edited by R. P. Behringer and J. T. 
Jenkins (Balkema, Rotterdam, 1997), p. 499. 



9 [m/s 2 ] 


1.62 


3.73 


9.81 


13.6 


25.1 


274 


en 


48.2 


41.4 


35.0 


33.3 


30.2 


18.4 


C/-R 


0.291 


0.274 


0.277 


0.269 


0.283 


0.260 



TABLE I. Angle of repose, G, in the drum middle and di- 
mensionless characteristic length, (/R, as function of gravity, 
g (simulation). 



G 



Figure Captions 

Figure |l| (a) Flat surface for low rotation speeds, (b) 
deformed surface for medium rotation speeds with 
two straight lines added as approximation and (c) 
fully developed S-shaped surface for higher rotation 
speeds. 

Figure ||| Experimentally measured starting (o) and 
stopping (*) angle and dynamic angle of repose (•) 
for mustard seeds. 

Figure |3| Dynamic angle of repose for black (•) and yel- 
low (o) mustard seeds. 

Figure |3] Dynamic angle of repose for small (•) and 
large (o) glass beads. 

Figure |H| Comparison of dynamic angle of repose for 
large mustard seeds taken from MM (o) and non- 
MRI (*) measurements. 

Figure |^ Comparison of dynamic angle of repose for 
large mustard seeds taken from MRI (o), numer- 
ical simulation (•) and the theory of Zik et al. |24|] 

(-)• 

Figure [7| Profile of the dynamic angle of repose along 
the rotation axis for 2.5 mm spheres: (•) simula- 
tion, ( — ) fit. 

Figure |§| Dynamic angle of repose as function of sphere 
diameter for Q = 20 rpm (simulation); (*) end cap, 
(•) drum middle, ( — ) arcus-tangent fit. 

Figure 9| Dimensionless range of boundary effect for 
spheres with different diameter; (*) R = 3.5 cm, (o) 
R = 7 cm, (•) R — 10.5 cm, dotted line shows value 
C = 0.28 R for radius-independent regime (simula- 
tion). 



Figure 10 Range of boundary effect measured in sphere 
diameters for spheres with different diameter; (*) 
R = 3.5 cm, (o) R = 7 cm, (•) R = 10.5 cm, (— ) 
exponential fit (simulation). 
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